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ABSTRACT 

Boxy and peanut-shaped bulges are seen in about half of edge-on disc galaxies. Comparisons 
of the photometry and major-axis gas and stellar kinematics of these bulges to simulations of 
bar formation and evolution indicate that they are bars viewed in projection. If the properties 
of boxy bulges can be entirely explained by assuming they are bars, then this may imply that 
their hosts are pure disc galaxies with no classical bulge. A handful of these bulges, including 
that of the Milky Way, have been observed to rotate cylindrically, i.e. with a mean stellar 
velocity independent of height above the disc. In order to assess whether such behaviour is 
ubiquitous in boxy bulges, and whether a pure disc interpretation is consistent with their stellar 
populations, we have analysed the stellar kinematics and populations of the boxy or peanut- 
shaped bulges in a sample of five edge-on galaxies. We placed slits along the major axis of 
each galaxy and at three offset but parallel positions to build up spatial coverage. The boxy 
bulge of NGC 3390 rotates perfectly cylindrically within the spatial extent and uncertainties 
of the data. This is consistent with the metallicity and a-element enhancement of the bulge, 
which are the same as in the disk. This galaxy is thus a pure disc galaxy. The boxy bulge of 
ESO 311-G012 also rotates very close to cylindrically. The boxy bulge of NGC 1381 is neither 
clearly cylindrically nor non-cylindrically rotating, but it has a negative vertical metallicity 
gradient and is a-enhanced with respect to its disc, suggesting a composite bulge comprised 
of a classical bulge and bar (and possibly a discy pseudobulge). The rotation of the peanut- 
shaped bulge of NGC 5746 is difficult to classify, but the peanut-shaped bulge of IC 4767 does 
not rotate cylindrically. Thus, even this relatively small sample is sufficient to demonstrate that 
boxy bulges display a range of rotational and population properties, indicating that they do not 
form a homogeneous class of object. 

Key words: galaxies: abundances — galaxies: kinematics and dynamics — galaxies: bulges 
— galaxies: stellar content 



1 INTRODUCTION 

The central regions of disc galaxies host bulges whose nature (or 
absence) must be explained by any model of galaxy formation and 
evolution. A bulge is a physical protrusion of stars above and below 
the plane of the disc and/or an excess of light above the inwards 
extrapolation of the radial exponential profile of the disc. There 
are three principal observational classes of bulges: classical bulges. 



* Based on observations collected at the European Organization for Astro- 
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0545, 65.N-0126 and 66.B-0073) 
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pseudobulges, and boxy or peanut-shaped bulges (see e.g. |Athanas-| 
|soula|2005| l. 

Classical bulges have dynamical and photometric properties 
and stellar populations that are similar to elliptical galaxies of the 
same mass ( Kormendy & Illingworth 1982 , Davies et al. 1983 , Em^ 



sellem et al. 2004, Kormendy & Kennicutt 2004^Falc6n-Barroso 



et all 1200 6: Thom as & Davies| |2006[ [MacArthur, Gonzalez & 



Courteau 



[2009j and are thought to be the end products of the 



same formation processes: either monolithic collapse (e.g. |Eggen, 
"Lynden-Bell, & Sandage 1962, Larson 1974, Carlberg 1984; Pipmo 
I& Matteucci,2004|l or multiple hierarchical mergers (e.g.,White & 



|Rees|1978[|Cole et al.|19"94{ [Thomas, Greggio & Bender|1999f 

Pseudobulges have embedded spiral structure, disc -like near- 
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Figure 1. Stellar line-of-sight rotation curves and velocity dispersion profiles for two non-cylindrically rotating edge-on galaxies from the SAURON sample 
(the fast rotator elliptical NGC 3377 and the classical bulge of the SO NGC 4570), constructed by extracting pseudoslits from the integral-field kinematics 
presented in jEmsellem et al.|j2004 |i. The stellar kinematics along the major axis are shown in black, and those from pseudoslits offset from, but parallel to, the 
major axis are shown in order of increasing distance from the major axis in red, green and blue. The offset distance z is shown. Data from one side of the major 
axis are shown as filled diamonds (joined by lines), and data from the other side as open diamonds (with error bars omitted for clarity). The vertical dotted lines 
at R = 14" for NGC 4570 indicate the effective radius R^, i.e. the approximate extent of the bulge (beyond which the figure is shaded grey). The effective radius 
of NGC 3377 is 38", i.e. beyond the radial range shown. Note that, within each bulge at a given R, v falls systematically with height above the major axis, i.e. 
the rotation pattern is not cylindrical. 



exponential light profiles, flattened axial ratios, and rotation- 
dominated dynamics (see [Kormendy & Kennicuttf 2004 and ref- 
erences therein). They are thought to be secularly rearranged disc 
material that was driven inwards by bars, ovals and possibly spi- 
ral arms. The observational constraints on the stellar populations of 
pseudobulges are less strong because studies have focussed on big- 
ger and brighter classical bulges, decomposing their extended star 
formation history is more complicated than in classical bulges and 
ellipticals, and dust is a significant problem for the late-type discs 
in which pseudobulges are more usually found. Nevertheless, the 
available population data are not strongly inconsistent with the sec- 
ular model (Peletier et al. 1999 , Gadotti & dos Anjos 2001 but see 
the pilot sample of ^Mac Arthur, Gonzalez & Courteau,2009| . 

An important third class of bulge, and the focus of this paper, 
is the boxy and peanut-shaped bulges, which are present in about 
half of edge-on disc galaxies jLtitticke, Dettmar & Pohlen'2000i 
Their major-axis stellar and gas kinematics |Kuijk en & Merrifield 
1995 , Merrifield & Kuijken 1999; Bureau & Freeman 1999 , Chung 
& Bureau.2004, Mendez-Abreu et al..2008j and hght profiles ^Bu-, 
reau et al. are consistent with simulations of the formation 

and buckling of galactic bars (Combes et al.'1990','Raha et al.'1 99l| 
[Athanassoula & Bureau 1999; Bureau & Athanassoula 2005] and 
point to them being bars viewed in projection. Peanut-shaped bulges 
are thought to be bars whose orientation is exactly perpendicular to 
the line-of-sight, while boxy bulges are seen in galaxies whose bar 
is at an intermediate angle. Bars oriented exactly parallel to the line- 



of-sight do not appear boxy or peanut-shaped. In that sense, boxy 
and peanut-shaped bulges are not the axisymmetric features nor- 
mally thought of as 'bulges'. Since the dominant model is that they 
are bars, which are not the end products of collapse or merger, but 
rather redistributed disc material, boxy bulges are also sometimes 
referred to as pseudobulges. However, for the purposes of clarity in 
the present paper, we avoid this, since boxy bulges are truly thick. 

In the absence of transformational mergers, bars, which are 
found in two-thirds of disc galaxies ^Eskridge et al.|2000[ |Whyte| 
|et aI.|20(j2{|Marinova & Jogee|2007^ , are expected to play a crucial 
and perhaps dominant role in disc galaxy evolution (e.g. |SelIwood 
& WiIkinson|p^93| [Kormendy & Kennicutt|[2004{ [Hopkins et al. 
20I0[ |. Boxy and peanut-shaped bulges therefore provide a unique 
perspective on a feature that is crucial for understanding galaxy evo- 
lution. By observing them above the plane of the disc, one can de- 
termine the intrinsic properties of bars in a way that is almost free of 
thin disc light pollution. Moreover, the observed dynamics of bars 
as a function of height above the plane is a crucial additional con- 
straint on models of the nature of bulges, including ruling out the 
presence of a dynamically hot classical bulge/merger remnant. For 
example, I Shen et aLlj2010 t modelled the line-of-sight mean veloc- 
ity and velocity dispersion of M giant stars observed by the BRAVA 
survey in the Milky Way bulge ( [Howard e t al. 2009). They placed 
an extremely low upper limit on the mass of any hot spheroid com- 
ponent, <7%; i.e. the Milky Way appears to be a pure disc galaxy 
(but see [Babusiaux et ar][2010[ l. Giant pure disc galaxies present 
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an acute challenge for cosmological simulations of galaxy forma- 
tion (e.g. [van den Bosch|2001[ [CmoUo et al.|2007[ |Dutton|2009[ 



IKormendy et al.|2010) . 

Perhaps the most unusual property of dynamical models of 
bars viewed in projection is that, within the bulge, the rotational 
velocity does not change with height (e.g. [Combes et al.|[T990[ 
[Athanassoula & Misiriotis|2002^ . That is to say, lines of constant 
velocity are parallel to each other and perpendicular to the major 
axis, such that (9v/i9|z| ~ 0, where v is the line-of-sight velocity 
and z is along the minor axis of the (edge-on) galaxy. This unusual 
behaviour was first observed in the boxy bulge of NGC 4565 by |Ko-| 
|rmendy & Illingworth| ( |1982t , who named it 'cylindrical rotation'. 
It has since been noted in the following eight galaxies, all of which 
have a bulges that is, more or less, boxy: NGC 128 ( |Jarvis|1990^, 
NGC 3069 ( |Shaw, Wilkinson & Carter 19931, NGC 1055 ( |Shaw 



1993), NGC 4442 I jBettoni & Galletta 1994j), NGC 7332 jPisher, 



niingworth & Franx|1994'{[Farc6n-Barroso et al.|2004^ , NGC 4220 
and NGC 4425 jFalcon-Barroso et al.||20 06), and the Milky Way 
( [Howard et al.|2009^ . In simulations, all bulges that rotate cylindri- 
cally are boxy or peanut-shaped, but not all boxy or peanut-shaped 
rotate cylindrically ( [Athanassoula & Misiriotis|2002| l. Cylindrical 
rotation is not present (or has not been reported) in simulated un- 
barred galaxies, or in observations of classical bulges, where rota- 
tional velocity falls systematically with height within the bulge (see 
Fig.[T]for two illus trative examples from the SAURON sample; [de] 
[Zeeuw et al.|2002[|Emsellem et al.|2"004l l. 

In our search of the literature, we found no observations of 
the absence of cylindrical rotation in a boxy bulge. One of the 
goals of the present targeted study of boxy bulges is to determine 
whether this is confirmation bias, or whether cylindrical rotation is 
truly ubiquitous in boxy bulges and therefore a requisite property 
of all realistic models of bars. The extent or absence of cylindrical 
rotation in boxy bulges would be of interest because it might con- 
strain the concentration of the dark haloes ^Athanassoula & Misiri-| 
[otis|2002| i. 

Our second goal is to examine the stellar populations of boxy 
bulges. The null hypothesis might be that, since bars are secularly 
rearranged disc material (rather than material accreted or formed 
during rapid mergers), the metallicities of boxy bulges should be 
similar to those of their host discs. If, however, the metallicities of 
boxy bulges differ from those of their discs, or if they have metallic- 
ity gradients, this is not necessarily a problem for the bar buckling 
scenario. As pointed out by [Freeman| ( [2008) , the stars that have been 
scattered furthest from the disc are the oldest stars, and therefore 
formed from the least metal-enriched fuel. The buckling process 
may therefore establish a negative minor-axis metallicity gradient. 
Indeed, such behaviour is observed in the Milky Way ([ Zoccali et al.| 
[2008^ and in NGC 4565, the archetypal boxy, cylindrically-rotating 
bulge ( [Proctor, Sansom & Reid|2000 1. The strength of the metallic- 
ity gradient established in this way must be a function of both the 
rate of enrichment in the disc and the timescale over which the bar 
buckles and heats the disc. The process must also establish a cor- 
responding positive age gradient. We know of no chemodynamical 
simulations of this process and, in particular, no prediction of the 
extent to which minor-axis gradients in boxy bulges are expected to 
differ from those in classical bulges. 

Moreover, we know of no prediction regarding the relative 
abundance of cc-elements in boxy bulges. The naive expectation, 
however, is that the [a/Fe] ratio of a boxy bulge should match that 
of its host disc. It is difficult to understand how secularly rearranged 
disc material could be segregated such that a boxy bulge of disc ori- 
gin is a-enhanced with respect to its disc. Enhanced [a/Fe] ratios 
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Figure 2. A'-band images of the sample galaxies with the slits overlaid. Con- 
tours are separated by 0.5 mag. Images taken from jBureau et al.]j2006^ . 



indicate that material was a-enriched by core collapse (type II) su- 
pemovae before type la supernovae Fe-enrichment, which implies 
a formation timescale < 1 Gyr. Bulges that are a-enhanced with re- 
spect to their disc therefore contains stars formed on much shorter 
timescales than the disc stars (and not just at earlier times). If ob- 
served in a galaxy, this would suggest that its bulge is at least par- 
tially comprised of material formed in a merger or rapid collapse, 
or material of external origin. We discuss some of the issues that 
complicate this naive expectation in Section|5] 

In this paper we present a focussed spectroscopic analysis of 
the boxy bulges of five edge-on disc galaxies. In Section |2] we de- 
scribe the sample, the data and their reduction. We present the de- 
rived stellar kinematics in Section |3] and the stellar populations in 
Section H] We summarise our results and suggest future work in 
SectionfS] 



2 OBSERVATIONS AND DATA REDUCTION 

Our sample consists of five galaxies whose bulges were classified 
as boxy or peanut-shaped by [Liitticke, Dettmar & Pohlen] ( [2000| l 
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Table 1. Sample properties and slit positions. 



Galaxy RC3 Bulge type Bulge extent Offset slit positions 

classification (arcsec) (arcsec) 



ESO311-G012 


so/a? 


Boxy 


15 


6.0 


9.0 


12 


IC 4767 


S0+- 


Peanut 


18 


3.5 


5.25 


7 


NGC 1381 


SAO 


Boxy 


10 


6.0 


9.0 


12 


NGC 3390 


Sb 


Boxy 


10 


5.0 


7.5 


10 


NGC 5746 


SAB(rs)b? 


Peanut 


20 


13.5 


20.25 


27 



Notes. Bulge extent is defined here as the radius of the maximum scale height of the disc, as measured from the vertical photometric parametrizations presented 
in Bureau et al. (in preparation). The offset slit positions are the distance each slit was offset from the major axis of the galaxy. 



using Digitized Sky Survey images, and whose major-axis stel- 
lar/gas kinematics and stellar structure have been homogeneously 
studie d ([Bureau & Freeman|1999[ [Chung & Bureau|2004[ [Bureauj 
[et al.[[2006p . All of the galaxies are edge-on or nearly edge-on. 
NGC 5746 is the furthest from edge-on in the sample. It is inclined 
at 85°. 5, measured by fitting an ellipse to the prominent ring visible 
in near-infrared Spitzer Infrared Array Camera images (PI: Gio- 
vanni Fazio). 

We observed the sample using the New Technology Tele- 
scope's ESO Multi-Mode Instrument (EMMI) in long-slit spec- 
troscopy mode. The observations were made in 2000 January, May 
and October (programmes 64.N-0545, 65.N-0126 and 66.B-0073). 
We positioned the slit along the major axis and then parallel to it at 
three offset positions. In this way we built up spatial coverage sim- 
ilar to that of a sparse, wide-field integral-field spectrograph. We 
present the sample properties and slit positions in TablefT] We show 
the slits superimposed on images of the sample in Fig. [2] 

We observed three other galaxy bulges (NGC 1055, NGC 1247 
and NGC 7123), two of which are not boxy, during the runs but 
we dropped them from the present analysis because the data are 
poor or incomplete, and the SAURON project and other work have 
superseded any contribution these observations could have made 
to our understanding of the dynamics and populations of classical 
bulges (e.g.jProctor & Sansom|2002|[EmseIlem et al.|2004[[Thomas| 
[& Davies|2006[[Falc6n-Barroso et al.|2006[ l. 

We reduced the data using standard techniques in IRAF (Image 
Reduction and Analysis Facility), yielding flat-fielded, wavelength- 
calibrated sky-subtracted two-dimensional spectra at four slit posi- 
tions for each galaxy. Exposure times ranged from 30 minutes for 
the major axis slits to up to 3 hours for the slits furthest from the 
disc. The reduced two-dimensional spectra were linearly binned in 
wavelength and cover the range 4830-5470 A. The spectral resolu- 
tion of the spectra is 1.0 A full-width half-maximum (equivalent to 
Oinst ~ 25 kms^' at 5150 A). The spatial axis has a pixel scale of 
0.9 arcsec/pixel. We determined the central bin (i.e. the spectrum for 
which the cylindrical radius R = 0) using a Sersic fit to the profile of 
each spectrum. This measurement is correct at the sub-arcsec level 
for the major axis slit, but is uncertain by < 5 arcsec away from the 
midplane, where the radial profile is noisy and weakly peaked. The 
central bin for each slit may therefore be a point that does not lie 
exactly on the minor axis. This uncertainty may introduce a small, 
systematic relative offset between the kinematic profiles derived for 
each slit. 

Spectrophotometric flux standards were not observed during 
the observing runs, so that we are unable to precisely flux calibrate 
our data. However, as described in Section |4] this does not signifi- 
cantly affect our ability to measure line strengths. 



3 STELLAR KINEMATICS 
3.1 Methods 

In order to determine the stellar kinematics, we first binned the data 
spatially to enhance the signal-to-noise ratio (S/N). To measure the 
line-of-sight mean velocity v, we typically binned to a 5/A' of about 
10 per A, but lowered this to 5 per A for a few of the most offset 
slits. The outermost bin in each slit often does not reach the desired 
S/N level, but we analysed it (and included in these figures) for 
completeness. To measure the velocity dispersion a, we rebinned 
the data to twice the 5/A', sacrificing spatial resolution. 

We extracted absorption line stellar kinematics from the 
binned spectra using penalized pixel fitting (pPXF; jCappellari &| 
[Emsellem 2004 ). All galaxies except NGC 1381 show clear signs 
of H/3 or [O III] emission in their major-axis spectrum, but the data 
in each bin are of insufficient S/N to constrain the emission line 
properties and subtract the emission in a physically motivated way 
using, e.g., GANDALF (Sarzi e t al.|2006^ . To derive the absorption 
line stellar kinematics, the appropriate spectral regions are there- 
fore simply masked[^We also masked a series of bad columns on 
the CCD at a rest wavelength of Ri5050 A. 

We used a subset of 88 stars from the MILES library of 985 ob- 
served stellar spectra as templates (Sanchez-Blazquez et al. 2006|. 
The MILES spectra have a spectral resolution of 2.3 A while our 
spectra have a resolution of 1.0 A so we first degraded our spec- 
tra to the MILES resolution. In principle, the penalized pixel fit- 
ting algorithm can robustly recover dispersions down to around 
half of this degraded spectral resolution, i.e. pa28kms^'. Using 
the Fourier Correlation Quotient method jBender|1990|l and a tem- 



plate star w ith a spectral resolution of ~32kms , Chung & Bu 



j2004 ^ measured dispersions > 50 km s ' along the major axis 



of all five galaxies in the present sample. Since the true dispersions 
of these galaxies is of order the spectral resolution of the MILES 
library, working at the MILES spectral resolution should not af- 
fect our results for the line-of-sight velocity dispersion. To verify 
this, we repeated the kinematic extraction for NGC 1381 working 
at the spectral resolution of our observations. To do this we used 
degraded template spectra from the ELODIE library of 1388 stars 
(whose spectra have a resolution of I2kms^', i.e. higher than that 
of our observations; [Prugniel & Soubiran|200I[ >. Within the uncer- 
tainties, the V and a radial profiles we obtained agree with those 
recovered using the MILES library, with no sign of any systematic 
disagreement. 

Our data are of insufficient S/N to constrain /13 (skewness) and 
/j4 (kurtosis) using a Gauss-Hermite parametrization of the line-of- 



' Note that we do treat emission properly when determining the absorption 
line strengths in Section|4] where we bin to much higher S/N. 
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Figure 3. Stellar line-of-sight mean velocity (top) and velocity dispersion (bottom) for the five sample galaxies. The stellar kinematics of the major axis are 
shown in black, and those of the offset slits are shown in order of increasing distance from the major axis in red, green and blue. See Table[T]and Fig.|2]for 
the slit positions. The vertical dotted lines indicate the approximate extent of the bulge (see Table^. The shaded regions at large radii therefore correspond 
roughly to the galaxy disc. The same systemic velocity has been subtracted from all points, determined from a pPXF fit to the total major axis spectrum. 
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Figure 3. — continued. 



sight velocity distribution (LOSVD). We therefore enforce a pure 
Gaussian LOSVD. Our conclusions regarding cylindrical rotation, 
which depend only on v, are not affected by this decision. 

We estimated the uncertainties on the derived kinematics by 
running pPXF on 100 random realisations of each spectrum. The 
realisations were constructed by adding wavelength-independent 
Gaussian noise assuming Poisson statistics. To these uncertainties 
we added uncertainties due to wavelength calibration, which we es- 
timate to be < 0.2 A, i.e. < 10 km s^ ' . These wavelength calibration 
uncertainties generally dominate the uncertainties in v, but are less 
important for <J, where the spectral noise dominates. 



3.2 Results 

The line-of-sight mean stellar velocity v and velocity dispersion CT 
calculated as described above are presented in Fig. [3] The kinemat- 
ics of the major axis slit are shown in black and those of the offset 
slits are shown in order of increasing distance from the major axis 
in red, green and blue. In each plot, we indicate the approximate 
extent of the bulge with dotted vertical lines. 

The boxy bulge of NGC 3390 is a clear example of cylindrical 
rotation. Within the uncertainties, it shows no systematic variation 
of rotational velocity as a function of height at all. 

The behaviour of the boxy bulge of ESO 311~G012 is close 
to cylindrical. The major-axis rotation is somewhat discrepant from 
that of the offset slits; the rotation curve along the plane of the disc 
is 'double-humped', while those of the offset slits are not. This does 
not mean that the boxy bulge is not rotating cylindrically. The strong 
double-hump feature is seen in simulations of barred galaxies (e.g. 



[Bureau & Athanassoula|2005) and is likely the dynamical signature 
of bar-driven inward transport of material (e.g. |Kormendy & Kenni-| 
|cutt|2004l l, which leads to the formation of a flattened round central 
excess of material, i.e. a pseudobulge. The distinguishing feature of 
cylindrical rotation is therefore that v does not vary systematically 
with height away from the major axis, which is true of ESO 311- 
G012. Taken together, we interpret the presence of a double hump 
in the rotation curve along the major axis and cylindrical rotation in 
the offset slits as consistent with simulations of the formation and 
dynamics of bars i Combes et al. 1990{ |Athanassoula & Misiriotis| 
|2002[|Bureau & Athanassoula.2005i) . 

The classification of the type of rotation in the boxy bulge of 
NGC 1381 is marginal. The velocity along the major axis is ex- 
tremely discrepant from that along the offset slits. The slits at 6" 
and 9" are corotating, but the slit at 12" is of extremely low S/N, 
and is binned to a correspondingly coarse spatial resolution. While 
this galaxy lacks the clear systematic decrease in v with increasing 
z seen in the examples in Fig.[T] it is not a clean case of cylindrical 
rotation either. 

The same is true of the peanut-shaped bulge of NGC 5746, 
but the picture is further complicated here by the rather irregular 
rotation curves within the bulge. This may be due to the prominent 
dust 'lane' in NGC 5746, which, because this galaxy is a few de- 
grees away from edge-on obscures a fraction of the galaxy. Strong 
absorption may cause the local mean line-of-sight velocity to differ 
significantly from the azimuthal velocity at the tangent point along 
the line-of-sight, and the extent of this difference may vary in a 
complicated way with distance along the minor axis. 

The peanut-shaped bulge of IC 4767 is not cylindrically ro- 
tating. At negative R, even accounting for the uncertainties and the 
irregularity of the rotation curve, v falls systematically with z- At 
positive R, the signature of non-cylindrical rotation is less clear. 

Cylindrical rotation is therefore not a ubiquitous feature of 
boxy and peanut-shaped bulges. There is also some evidence that 
it is less strong or even absent in peanut-shaped bulges compared 
to boxy bulges. It is worth considering whether this might have a 
dynamical origin. In simulations, bars in edge-on galaxies appear 
peanut-shaped when viewed with the bar perpendicular to the line- 
of-sight, and boxy when the bar is at an intermediate angle (e.g. 
[Combes et al.|1990 l. End-on bars appear neither boxy nor peanut- 
shaped. The velocity fields of these simulated galaxies vary with the 
orientation of the bar in the sense that the closer the bar is to end-on 
(i.e. the less strongly peanut-shaped it appears), the more closely 
packed are the lines of constant velocity (e.g. figure 6 of | Combes| 
[etal.[ 1990 ' and figure 12 of Athanassoula & Misiriotis[2002[ l. There 
is, however, no clear sign of a systemic strengthening or weaken- 
ing of cylindrical rotation as the viewing angle varies in these pub- 
lished studies. Rather, the strength and extent of cylindrical rotation 
in the simulations of | Athanassoula & Misiriotis|j2002 ) is a function 
of dark halo concentration in the sense that less concentrated dark 
haloes host boxy bulges in which cylindrical rotation is absent or 
less extended. 

On a different note, in ESO 31 1~G012, NGC 1381, NGC 3390 
and NGC 5746 (i.e. all the galaxies with boxy or peanut-shaped 
bulges except IC 4767), there is a general downwards trend in ve- 
locity dispersion with increasing height above the disc. This is most 
clear in ESO 3 1 1-G012 and NGC 5746. This behaviour is also seen 
in classical bulges (Fig.[T|l, but the details of the variation of the ve- 
locity dispersion with height were crucial constraints in the Milky 
Way models of [Howard et al.[j2009) and [Shen et aEjpOTOl l. These 
authors placed an upper limit on the possible contribution of any dy- 
namically hot spheroid component by comparing photometric and 
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Figure 4. Example GANDALF fit for one of the two disc spectra of NGC 3390. We show the observed spectrum (degraded to the MILES resokition) in black 
and the GANDALF best fit in red. The green points are the difference between the best fit and the observed spectrum, shifted by an arbitrary positive offset so 
they appear on the scale. The best fit is comprised of the H/3 and [O III] emission lines (shown in blue) and a combination of MILES template stars that are 
convolved by a LOSVD (purple dashed line) and multiplied by a 6th order polynomial to account for the relative continuum differences between our setup and 
the MILES library. The anomalous feature at «4920 A is due to a problem with the CCD, and is masked for all analyses. 



kinematic observations to dynamical models. Their analysis implies 
that the Milky Way has, at most, an extremely low-mass classical 
bulge, i.e. it is a pure disc galaxy (but see |Babusiaux et al.|2010^ . 
Without a comparison to dynamical models custom-made for each 
galaxy, which is beyond the scope of the present work, it is not 
possible here to draw such a strong conclusion on the basis of our 
velocity and dispersion profiles. 



4 LINE STRENGTHS AND STELLAR POPULATIONS 
4.1 Methods 

The spatial variations of stellar population properties offer addi- 
tional constraints on the nature and evolution of bulges. We there- 
fore measured the strengths of the absorption lines present in our 
data (Hj3, Fe5015, Mgb, Fe5270, Fe5335 and Fe5406) in the 
Lick/IDS system jBurstein et a l. 1984 : Worthey|1994[[lYager et al.| 
1 1998) and compared these to single stellar population (SSP) mod- 
els, yielding SSP-equivalent luminosity-weighted ages, metallici- 
ties [Z/H] and a-element enhancement [a/Fe]. 

To make these measurements, we require high S/N to reliably 
subtract nebular emission and thus accurately measure the absorp- 
tion line properties. We therefore binned the major-axis slit into 
three large bins: a single bin for the bulge and two disc bins (one 
from each side of the galaxy). We collapsed each of the three offset 
slits into a single bin dominated by bulge light (rejecting the very 
small amount of thick disc light at larger galactocentric radii). This 
yields one bulge-dominated measurement from the centre of the 
bulge, two disc-dominated measurements along the major axis and 
three minor-axis bulge measurements at increasing galactic height. 
With data at four positions along the z-axis, it is not really pos- 
sible to quantify vertical gradients (although as we discuss in the 
introduction, these are of interest). However, we can test our naive 
picture that a boxy bulge should have [a/Fe] similar to that of its 
host disc. 

Of the absorption lines in our data, only Hj8 provides a strong 
constraint on stellar age that is weakly dependent on metallicity 
( [Worthey]|1994| ), but this absorption feature is unfortunately of- 



ten filled in by emission. To compare the population properties of 
galaxy discs and bulges, we therefore restrict our analysis to those 
galaxies whose major-axis disc spectra are either free of emission 
or of sufficient S/N to allow reliable subtraction of the Hj3 emis- 
sion using GANDALF (Sarzi et al. 20061. This is only true for 
NGC 1381, which has no emission, and NGC 3390. The analysis 
that follows is therefore restricted to those two galaxies only. 

We show an example GANDALF fit in Fig.|4] The best fitting 
GANDALF solution is comprised of the Hj3 and [O III] emission 
lines and a combination of flux-calibrated MILES template stars 
that are convolved by a LOSVD and multiplied by a 6th order poly- 
nomial to account for the relative continuum differences between 
our setup and the MILES library. To produce an emission-cleaned 
spectrum that can be used to measure absorption line properties and 
that can be compared to models, we subtract the GANDALF emis- 
sion lines from the observed spectrum and divide it by the 6th order 
polynomial, achieving a rough relative flux calibration. In practice, 
the polynomial is extremely slowly varying and close to first order 
over the bandpasses of the Lick indices, so this affects the measured 
line strengths at a level significantly below the uncertainties due to 
the noise. 

It is then possible to measure the strengths of the lines in the 
Lick/IDS system, which we do using the method of |Kuntschner| 
( [2000[ >. But before these measurements are made, the spectra must 
be convolved with a wavelength-dependent Gaussian to account for 
the difference between the spectral resolution of our setup and that 
of the Lick/IDS system. The resulting indices are further corrected 
for (1) the broadening effects of the galaxies' velocity dispersions 
(which generally reduce the strength of absorption features) and (2) 
the small systematic offsets caused by, e.g., the remaining differ- 
ences between the continuum shapes our spectra and the Lick/IDS 
spectra. We apply the offsets from the MILES system given in the 
appendix of |Paudel et al.|(2010| l. This offset of course only affects 
the absolute value of the indices, not their spatial gradients. 

We use the line strengths to measure stellar population prop- 
erties using the models of [Thomas, Maraston & Bender| p003) . 
These models incorporate non-solar [a/Fe]. The best-fitting SSP- 
equivalent model for each spectrum is found by exploring an inter- 
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NGC 1381 



NGC 3390 



polated grid of model Lick indices in age, [Z/H] and [a/Fe] (see, 
e.g., |Proctor, Sansom & Reid|2000[|Kuntschner et al.|2010[ l. We ex- 
clude Fe5015 from this analysis. It is still sometimes contaminated 
by [O III] emission, it is close to the bad pixels on the CCD, and it 
covers an extremely wide bandpass so is subject to relatively large 
continuum systematics and Lick offsets. This line provides very lit- 
tle additional constraints on the optimal SSP over those offered by 
the other three Fe lines. 



4.2 Results 

We show the derived Lick indices and corresponding SSP- 
equivalent luminosity-weighted age, [Z/H] and [a/Fe] as a func- 
tion of height above the galactic plane in Fig. [5] There are two addi- 
tional points on the major axis (z = 0), which are the measurements 
from the two disc bins. 

In both galaxies, there is no sign of a systematic trend in the 
age-sensitive Hj3 line as a function of z, or in the derived age. 

In NGC 1381 there are broadly decreasing trends with increas- 
ing height in the metallicity-sensitive Fe and Mg lines. These are 
qualitatively similar to the minor-axis trends seen in the SAURON 
classical bulges shown in Fig. [T] (see |Kuntschner et al.||2006| for 
the corresponding line strength maps). The SSP-equivalent stellar 
population properties derived from these line strengths suggest that 
the bulge is comprised of an old stellar population, > lOGyr with 
no systematic trend in age as a function of z, but there is a clear 
decreasing trend in [Z/H]. We also note that [a/Fe] in the bulge at 
large z is greater than that in the disc bins (and indeed the central bin 
of the major axis). The simplest interpretation of these results is that 
NGC 1381 is comprised of a disc and a bulge whose stars formed 
rapidly. This explains the general trend in [a/Fe] as a function of 
height, as disc light (with its lower [a/Fe]) contributes less and less 
to the integrated spectrum. In summary, the bulge of NGC 1381 is 
chemically similar to that of classical bulges and elliptical galaxies. 

Given this classical bulge behaviour, how do we explain the 
boxy appearance of the bulge? We speculate that this is a compos- 
ite bulge: its chemical properties are explained by the presence of a 
classical bulge, but its appearance suggests the simultaneous pres- 
ence of a bar (which appears boxy in projection). This is consis- 
tent with the stellar rotation, neither cylindrical nor strongly non- 
cylindrical. The double-hump of the rotation curve even hints at the 
presence of a small discy pseudobulge. This galaxy thus seems to 
have all three kinds of bulges. 

The situation in NGC 3390 is much simpler. Within the (not 
insignificant) uncertainties, its bulge has no minor-axis metallicity 
gradient, and the [a/Fe] measurements in the bulge bins at large 
z are entirely consistent with those of the disc. That is to say, this 
galaxy's bulge is made of material similar to that in its disc, and the 
process that has scattered it to large z has not established a system- 
atic metallicity gradient. This is consistent with the kinematics of 
this bulge, which is cylindrically rotating to great precision in all 
the slits (including the major axis). The simplest interpretation of 
NGC 3390 is therefore that it is a pure disc galaxy. 



5 DISCUSSION 

We analysed the kinematics and stellar populations of the boxy or 
peanut-shaped bulges of a sample of five edge-on disc galaxies. One 
boxy bulge (NGC 3990) rotates perfectly cylindrically, while an- 
other (ESO 31 1-G012) is very close to cylindrical. It is not possible 
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Figure 5. Lick indices and luminosity-weighted SSP-equivalent ages, 
metallicities and [a/Fe] ratios as a function of the perpendicular distance 
to the major axis for NGC 1381 and NGC 3390. Measurements are made at 
the center of the galaxy (black squai'e at z = 0), two disc bins along the major 
axis (blue diamonds) and three points at increasing distance along the minor 
axis (black squares at z ^ 0). The disc bins correspond to spectra at z = 0, but 
we added a spatial offset of ±1" for clarity in these plots. The z = 12" bulge 
spectrum in NGC 1381 has some line strengths that are consistent with an 
equivalent width of zero, and the uncertainties in [a/Fe] for this bin extend 
over the full model grid. The position is shown merely for completeness. 
The Fe5406 feature is redshifted beyond the observed wavelength range for 
NGC 3390 and could therefore not be measured. 
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to cleanly classify the rotation in the boxy bulge of NGC 1381 or 
the peanut-shaped bulge of NGC 5746. The strongly peanut-shaped 
bulge of IC 4767 is in normal rotation, the velocity falling systemat- 
ically with height. Cylindrical rotation is therefore not a ubiquitous 
behaviour in bulges that appear boxy or peanut- shaped. 

Our sample is small and this work is exploratory, but we also 
tentatively note that cylindrical rotation may be less strong or even 
absent in peanut-shaped bulges, compared to boxy bulges. While 
this is not a prediction of simulations, models are able to qualita- 
tively reproduce a case like IC 4767 by reducing the central con- 
centration of the dark halo ( [Athanassoula & Misiriotis|2002) . 

Low 5/^ and the presence of emission are less of a prob- 
lem when determining absorption line stellar kinematics than when 
measuring absorption line strengths. We were therefore only able to 
measure minor-axis trends in Lick indices and derive luminosity- 
weighted SSP-equivalent age, [Z/H] and [a/Fe] for two of the 
five galaxies. Of these, NGC 3390 exhibits behaviour consistent 
with being a pure disc: there is no significant difference in [Z/H] 
and [a/Fe] between the disc and the boxy bulge at large galactic 
heights, suggesting that the boxy bulge formed from redistributed 
disc material. NGC 1381, however, shows a clear negative metallic- 
ity gradient as a function of z, while [a/Fe] in the bulge at large z is 
significantly enhanced compared to that in the disc. This behaviour 
is consistent with a simple disc + classical bulge model (e.g. |Nor^ 
|ris. Sharpies & Kuntschner|2006| for the edge-on SO NGC 31 15). 
Combined with the marginal cylindrical rotation and the boxy shape 
of the bulge, this suggests that the centre of NGC 1381 hosts a com- 
posite bulge comprised of a classical bulge, a discy pseudobulge and 
a thick bar. 

The naive expectation that we described in our introduction, 
that a boxy bulge should have the same [a/Fe] as its host disk if 
it formed through the bar buckling scenario and is part of a pure 
disk galaxy, is one part of our argument that NGC 3390 may be a 
pure disk galaxy. This picture is complicated by at least two issues. 
(1) The disk could have formed rapidly and then buckled promptly, 
forming an a-enhanced thick bar, i.e. a boxy bulge, and (2) In a sys- 
tem made up of stars in which [a/Fe] is anticorrelated with [Z/H] 
(e.g. as seen in young stars in the Galaxy, [Edvardsson et al.|1993^ , 
a negative spatial metallicity gradient would necessarily coincide 
with a positive spatial [a/Fe] gradient. 



Regarding point (1), we note that the 'wild disc' scenario 
(in which early star formation in the progenitors of disk galaxies 



takes place on short timescales in massive clumps, e.g. 


Noguchi 


1999: Elmegreen, Elmegreen & Hirst|2004| |Boumaud, E 


megreen 


& Elmegreen 2007 1 probably cannot explain how a bar could be 



a-enhanced with respect to its host disc. In dynamical simulations, 
these clumps rapidly migrate inward through dynamical friction and 
combine to form an axisymmetric feature that would probably look 
like an a-enhanced classical bulge (both in terms of populations 
and dynamics). Chemodynamical simulations of bar formation are 
clearly necessary to clarify the issue. 

Such studies should also address the more general questions 
raised by this work. (1) Is the double hump we see in the rotation 
curve along the major axis of two of the boxy and peanut-shaped 
bulges definitely associated with a flattened discy pseudobulge? 
How far does its influence extend above the plane? (2) How does 
the viewing angle of the bar affect the observed kinematics? For ex- 
ample, do we expect weaker cylindrical rotation in peanut-shaped 
bulges? This question was addressed by Atha nassoula & Misiriotis| 
( |2002[ >, but quantitative (rather than qualitative) comparisons to ob- 
servations are necessary. (3) Do the minor-axis metallicity and age 
gradients of simulated bars differ systematically from those of clas- 



sical bulges? (4) Can simulated bars form in a way that explains 
the presence of a positive minor-axis [a/Fe] gradient, as seen in 
NGC 1381, without the additional presence of a classical bulge? 

The a-enhanced bulge of NGC 1381 suggests that not all boxy 
bulges are found in pure disc galaxies (although some, including 
the Milky Way, may be). Since the prevalence of giant pure disc 
galaxies is an acute problem for simulations of galaxy formation in 
a cosmological context, it is clear that a combined analysis of the 
stellar dynamics and populations of the central regions of galaxies 
is an extremely powerful tool. 

The sample we have studied in the present work is unusual 
because it is dominated by boxy and peanut-shaped bulges. Such 
galaxies are unique laboratories to investigate the nature of the bars 
that drive secular evolution. Moreover, they raise the possibility that 
a demographic ally significant fraction of giant disk galaxies, mak- 
ing up about half the population, is composed of pure disc galaxies. 
We have shown that a property as simple as the degree of cylindri- 
cal rotation varies significantly among boxy/peanut-shaped bulges, 
as do the properties of their stellar populations. The driver of this 
variation is unknown, but may be a property as fundamental as dark 
halo concentration. Clearly, if we want to further our understanding 
of the formation of boxy and peanut-shaped bulges, we need data in 
larger quantities (to look for correlations with other galaxy proper- 
ties such as environment and mass distribution) and at higher 5/A' 
(to more accurately measure stellar population properties). Two- 
dimensional data from integral-field units would also be beneficial, 
allowing the dynamics to be modelled in detail. 
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